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' We investigate the symmetry properties for Baym's ^-derivable schemes. We show that 

in general the solutions of the dynamical equations of motion, derived from approximations 
' of the ^-functional, do not fulfill the Ward-Takahashi identities of the symmetry of the un- 

derlying classical action, although the conservation laws for the expectation values of the 
corresponding Noether currents are fulfilled exactly for the approximation. Further we prove 
that one can define an effective action functional in terms of the self-consistent propagators 
which is invariant under the operation of the same symmetry group representation as the 
classical action. The requirements for this theorem to hold true are the same as for pertur- 
qq ' bative approximations: The symmetry has to be realized linearly on the fields and it must 

be free of anomalies, i.e., there should exist a symmetry conserving regularization scheme. 
In addition, if the theory is renormalizable in Dyson's narrow sense, it can be renormalizcd 
. with counter terms which do not violate the symmetry. 

<N 
O 



PACS numbers: ll.10.-z, ll.10.Gh, ll.10.Wx 



P-t. I. INTRODUCTION 

Symmetry principles are of fundamental importance in physics. For instance the standard model 
of elementary particle physics rests on the principle of local gauge symmetries. These symmetries 
do not only guarantee the observed conservation laws but also ensure the physical consistency of 
the quantum field theory. Especially no unphysical degrees of freedom for the vector particles 
appear falsely in the particle spectrum. This ensures the causality and the positive defmiteness of 
probabilities (in the vacuum case especially the unitarity of the S-Matrix). 

Though fulfilled for the exact theory such symmetry concepts are commonly violated in most of 
the approximation schemes, in particular if partial resummations are performed. In this paper we 
address the question of global symmetries for Dyson resummation schemes based on the two-particle 
irreducible (2PI) action formalism |l|, [2|, |j also known as ^-derivable approximations Jl. Here $ 
denotes the 2PI part of the generating functional T. As Baym could show such approximations 
have the property that expectation values of the Noether currents are exactly conserved and the 
solutions are thermodynamically consistent. These properties even survive a consistent gradient 
approximation of the Kadanoff-Baym equations leading to generalized quantum kinetic equations 
[]|] which dynamically treat the full spectral functions of the particles and therefore permit a 
consistent transport treatment of resonances. 

In the first paper of this series H, in the following referred to as [I], we could verify an other 
important property of this type of approximations, namely that of renormalizability, provided the 
original theory is renormalizable in the standard sense. Using the path integral techniques || 
within the Schwinger-Keldysh real time formalism (for details and the notation conventions used 
here see [I]) closed equations for the renormalized self-consistent quantities with a counter term 
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structure solely denned at the vacuum level could be established. First examples which include the 
fully self-consistent approximation at second order in the coupling, i.e., including the tadpole and 
the sunset self-energy diagrams, were given in the second part of this series 0, referred to as [II]. 

While for the solution to ^-derivable approximations the expectation values of Noether currents 
are exactly conserved, this is no longer guaranteed at the correlator level or for higher order 
vertex functions. Thus in general the Ward-Takahashi identities (and thus also the symmetries 
of the underlying classical action) are violated already for the self-consistent self-energy. This 
was seen by Baym and Grinstein investigating various approximations for the 0(iV)-model [||. 
Considerations related to symmetries in approximation schemes were conducted already very early, 
e.g., in investigations of superfluid Helium within non-relativistic many- body theory g . At that 
time the notion of gap-less approximations was coined in the context of Nambu-Goldstone modes 
in spontaneously broken phases. In some recent attempts the same phenomena were investigated in 
pi comparing various regularization and renormalization techniques, while in [11] the equations 



of motion were modified by introducing mean-field dependent effective vertex functions in the 
approximations at the gap-equation level. 

In this paper we like to investigate the symmetry aspects of approximation schemes at a quite 
general level fully within the 2PI functional approach. In Sect. II we derive generalized Ward- 
Takahashi identities for the 2PI action functional. This is done in close analogy to the corresponding 
investigations for the usual perturbative quantum action functional (see for instance [12]). Thereby 



we prove the following properties of the 2PI action functional formalism: Provided the symmetry 
is realized as a linear representation on the fields and is free of anomalies the considerations show 
the following: 

(1) If the 2PI action is truncated in a systematic way (loop expansion, 1/N expansion) it fulfills 
the same symmetries as the classical action if both, the classical fields and the Green's 
functions, are simultaneously transformed as tensors of the symmetry representation of 1 st 
and 2 nd rank, respectively. 

(2) For the solutions of the corresponding self-consistent Dyson equation: 

(a) expectation values of Noether currents are conserved. 

(b) The scheme is void of double counting and fulfills detailed balance relations and is thus 
thermodynamically consistent. 

(c) However, Ward-Takahashi identities at the correlator level and for higher order vertex 
functions can be violated. 

The latter problem is expected for such partial resummation schemes since from the point of 
view of the exact theory, which fulfills the symmetry, this expansion is of course incomplete. 
Diagrammatically it is easy to see that in general ^-functional based partial Dyson resummations 
already violate the crossing symmetry of vertex functions in orders of the expansion parameter 
beyond that chosen for the truncation of the ^-functional. 



As an important result of this paper it will be shown in Sect. Ill how to overcome these prob- 
lems. The strategy is that for any ^-derivable approximation based on a truncated 2PI action 
T[(p, G] one always can construct a non-perturbative approximation for the effective quantum ac- 
tion on top of the self-consistent solution which indeed recovers the symmetry in the above cited 
sense. This effective action, called T[ip] is symmetric. It is expressed in terms of self-consistent 
propagators G[np] at a given mean field and provides a non-perturbative approximation to the usual 
1PI quantum-action functional. This action was already considered in ref. Q where the identity 
of vertex functions derived from T and T has been shown for the exact case. Here we concentrate 



3 



on the consequences of truncation schemes, where both pictures are no longer equivalent. Rather 
in terms of perturbative diagrams, the new action indeed supplements that minimal set of extra 
diagrams needed to recover the symmetry for the restricted set of diagrams resummed by the <J>- 
derivable Dyson scheme. The extra terms are encoded in a special Bethe-Salpeter equation and 
higher vertex equations which precisely cope with the chosen ^-derivable scheme. Formally the 
functional derivatives of this effective action with respect to the background field, taken at the 
stationary point, define non-perturbative approximations to the self-energies and proper vertex 
functions which fulfill the usual Ward-Takahashi identities. We call the so generated vertex func- 
tions external, as they do not take part in the self-consistent scheme but are rather constructed a 
posteriori, once the self-consistent Dyson solutions are given. In this way symmetry-conserving non- 
perturbative approximations to these quantities are provided. It is shown that this also holds true 
for the renormalized physical quantities. Since the external vertex functions are defined as usual 
by multiple functional derivatives of the effective action all symmetries inherent in this functional 
like crossing symmetry are guaranteed for the solutions. We shall concentrate on the properties of 
the self-energy, which, if derived from T, respects the symmetries (e.g., posses Nambu-Goldstone 
modes). This general scheme, however, permits to repair the symmetry of any other higher order 
vertex functions such as the correlators of Noether currents {j tl (x)j l/ (y)) , which then are conserved 
in contrast to the ones given by the underlying self-consistent Dyson resummation. 

In Sect. IV we discuss different approximation levels for the 2PI action applied to the 0(7V) 
linear cr-model. Besides the Hartree approximation (truncating the functional at the order A) we 
also discuss the next approximation level, i.e., up to second order in the coupling. As the most 
simple example for the recovery of the 0(A r )-symmetry violated by the Hartree approximation in 
the broken Nambu-Goldstone phase we show that the well known Random-Phase approximation 
(RPA), see, for instance, [p~3|| , is derivable from a non-perturbative effective action of the here 
defined kind. Some explicit numerical results are given for the chiral linear <r-model. This clarifies 
the question of symmetry violations of the Hartree approximation shown in from the general 
point of view elaborated in this paper, which specifies the route to symmetry preserving vertex 
functions for any kind of ^-derivable approximation. In the context of this series it is shown in 
sect. Ill <H that a straight forward generalization of the renormalization procedure presented in 
the first part of this series [I] also leads to closed renormalized equations for these vertex functions 
at finite temperature which counter-term structures are solely defined on the vacuum level. The 
paper is closed with conclusion and outlook. 

The appendices contain a formal derivation of the expression for the ^-functional which uses 
the path-integral formalism developed in || and a short summary about Noether's theorem in the 
classical context which is needed in Sect. II to formulate the generalized Ward-Takahashi identities. 



II. SYMMETRY PROPERTIES OF THE EXACT 2PI ACTION FUNCTIONAL 

The 2PI action functional 0, §, £§ can be introduced with help of the well known path integral 
methods ||. Besides the standard local source J one also supplements a bilocal source B to 
the partition sum Z = Z[J,B]. The 2PI action functional T[(p, G] is then defined as the double 
Legendre transformation of W[J, B] = — \h\n Z[J, B] with respect to both, J and B. Details and 
a slightly different derivation than in |J are given in Appendix [A|. A formal loop expansion of the 
path integral yields 

r[^,G]=5M + |Trln(Gr 2 1 /M 2 ) + | f d(12) &£{Gn - + G]. (1) 

J c 

Here we have used the notation J c d(l) f± for the integration over cf-dimensional space time in the 
sense of dimensional regularization and appropriate sums for the internal field components (i.e., 
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"charge space" indices). Within the real-time formalism adapted to the equilibrium case the time 
integration runs along a modified Schwinger-Keldysh contour 1 which includes an imaginary branch 
of length i/3. In (p]) S[ip] denotes the classical action as a functional of the mean field tp, while 
the terms proportional to h account for the quantum fluctuations at the one loop level. The free 
propagator in the presence of the mean field <p is denoted by £): 



(5T 1 ) 



12 



(2) 



Finally $ accounts for all higher order terms. It is diagrammatically defined as the sum of all 
two-particle irreducible (2PI) closed diagrams beyond the one-loop level which can be built with 
point-vertices defined by the interaction part of S[6 + <p], i.e., the part with at least three ^-fields 
expanded in the same way as the Feynman diagrams of perturbation theory. Thereby lines stand 
for full propagators G rather than free propagators. We shall give an analytical definition of the $ 
functional which can be used for practical calculations of approximations like the coupling constant 
or loop expansion for the ^-functional in appendix |A|. 

The equations of motion for the mean field ip and the Green's function are given by the station- 
arity of T, i.e., at vanishing auxiliary sources J = B = as 



8<p ' 5G 

Using Eq. (|]) the second equation is seen as the Dyson equation for the full propagator 



(^~ 1 )i2 - (0 



12 



2i 



5G 



-42- 



12 



(3) 



(4) 



Thus (|) gives a closed self-consistent set of equations of motion for the mean fields and the self- 
energies in terms of the exact Green's function G. The 2PI property of $ avoids double counting 
in the sense that the diagrams of the self-energy with lines denoting exact Green's functions do 
not contain any self-energy insertion in any of its lines, i.e., it generates one-particle irreducible 
(1PI) skeleton diagrams for the self-energy by variation with respect to G. This is immediately 
clear from the fact that the derivative of $ with respect to G diagrammatically implies to open 
anyone of its lines in the diagrams building <I> and taking the sum over the so obtained diagrams 
with two truncated external points. 

For the discussion of the symmetry properties of the above defined functionals we take the 
0(./V)-symmetric (/> 4 -theory as an example. The generalization to other models and symmetries 
with more complicated field configurations is straight forward. 

We use the fact that the path integral measure is invariant under a field translation 6' = 6 + 56. 
In first order of 86 we find 



B6 



x cxp 



d(l) ^T + Jii 



% + J ^{12) B jhk2 6{86 k 2 



i 



iS[6) + i / d(l) J l3 6{ + d(12) B jltk2 6{6 k 2 



(5) 



besides the real-time parts along the real time axis from an initial time U to a final time tf and back to U this 
contour also includes a branch parallel to the imaginary time axis from U to fcj — i/3 (/3 = 1 /T denotes the inverse 
temperature of the system). Finally we let the times U and tf go to — oo and +oo, respectively. All considerations 
can be extended to more general off-equilibrium initial statistical operators with the qualification that in this case 
it does not make sense to take U — > — oo 
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Here subscripts and superscripts j, k denote the field components {(ft) = (ft, and Einstein's summa- 
tion convention is implied. For a local 0(iV)-transformation 



b{ = i5 X i(T a y r <ft{ 



(6) 



we obtain 
= 



hi 



d(i) 



5S 



+ J y I i(r a y 8x1 + 2 / d(2) B jl>k2 i{ry r 448 X l 



x exp 



;' 1 



iS[cft] +i / d(l') J Vj ,4? v + ^ / d(l'2 r ) B Vj ,p k ,<lP v <j%, 



(7) 



Since the classical action functional is invariant under global O(N) transformations (^), c.f. Ap- 
pendix [A], we read off from (B8) that @ can be expressed in the form 



d(l) 5x1 d^/ 



i5J 



+ 2 /d(i2) j^^y 



(8) 



Since <$Xi is an arbitrary function this can be brought to the local expression 



= d,jl 



a 1 1 



iSJ 



Z[J,B] + Jlj i(ry/-^- + 2 [ d(2) B ljj2k Kr a ) j f ^f^ 
3 SiJtf J c 3 5iB lfj2k 



(9) 



For the solution of the equations of motion we have J = and B = 0. Thus the expectation value 
of the current is conserved as we read off from (||): 



Clf.1 



_8_ 

i5J 



Z[J,B] = Z[J,B]d ll (ff) = 0. 



(10) 



Here the expectation value has to be read as the quantum statistical expectation value for the local 
current operator if interpreted within the operator formalism of quantum field theory. Here and 
in the following we write operators as bold face upright symbols. 

For sake of completeness we write down the local Ward-Takahashi identities (WTIs) for the 
current. From (|B11| ) we find for the 0(iV)-Noether current the expression 



fiM = -i(T a )n>4>{d»(i>{ • 
Using this in (S) yields the local WTIs for the functional W = —ilnZ: 



□ 



(2). 



5W 



5Bi k 



lk,2k' 



1 = 2 



1/ 



5W 



3 SB ritk2 



(11) 



0. (12) 



Finally we also obtain an expression for the r-functional by using the relations ( |A5| ) and ( |A 15 ): 



-^□y>f + 2i 



□ ( 2 )G 



lfc,2fc' 



1=2 



sr 

6(p 



^fVi + 2 / d(2) j^—r^G 
Pi Jc oGji, k2 JJ 



{T a )kk> 
lj',2k- 



(13) 



Again it follows immediately that for the solutions of the equations of motion (||) the expectation 
value of the 0(iV)-Noether current is conserved. Taking further derivatives of (|9|) yields WTIs for 
higher order Green's functions of the Noether current, which we do not need in the further line of 
arguments. 

For the further symmetry analysis we like to concentrate on properties of the 2PI functional (jl|). 
We thus go back to (0) and set 5x % = const. Then the first term in the curly bracket is a complete 
divergence and thus the integral vanishes and due to the linear independence of the constants Sx 1 
we obtain the global form of the generalized Ward- Takahashi identities. 

Especially this proves that the generating functional Z is invariant under global O(N) transfor- 
mations when the local source J\j is transformed in the contragredient way of the fields and the 
bilocal sources B\^ 2k as a covariant tensor of 2 nd rank. Since (§) is of first order in the derivatives 
of Z with respect to the sources the same holds true for W: 



1) JiAry / W s [ j: B] + /. d(12) B ljl!2k (ry' j + B lj>2kl (r a f k 



if 



5W[J,B] 
SBji jk2 



and from (A5) we find (again for global transformations): 



y c d(l) -^-r ^ 1 + y c d(12) _g- [ (T V f G 12 + (r ) k ,G l2 



(14) 



(15) 



This result can be derived also directly from by taking the integral J c d(l) ... on both sides 
of the equation. 



III. SYMMETRIES OF 2PI ^-DERIVABLE APPROXIMATIONS 

So far the Ward-identities were derived for the exact Junctionals, Green's functions and mean 
fields. In particular it is important to realize that the identity 

, k 5 2 W\J,B] 5W\J,B] , j. r , 

\G\ 2 := — i— = 2— ( / 3 iV 9 2 ^ or the exact functional (16) 

6J\j5J 2k oB\j y2k 

was proven using the underlying path integral definition of the exact functionals. We now step 
towards the ^-derivable approximations. They are defined by a truncation of the auxiliary func- 
tional $[ip, G] in the definition (||) of the generating functional T[ip, G] while keeping the variational 
properties (||) which define the equations of motion for mean field and propagator 

£r a p prx [y,G] ^ Q (5r apprx [v3,G] ^ Q 
5ip 5G 

The latter equation of motion defines the self-energy in terms of the self-consistent propagator by 

S 12 := (St^n - (G-% = 2i ^%^l. (18) 



5G 



12 



Thereby $[93, G] can be truncated according to various schemes, like expansion in powers of A, in 
loop order or in powers of 1/N (see |l4|]), excluding the internal structure of the Green's functions 
from the counting. The sole requirement for the thus constructed <3? a pprx is that it remains invariant 
under the symmetry transformations as explained in appendix [A], i.e., that the generalized WTIs 
(|l~5| ) hold true also for the approximation. In our case of the linear 0(A r )-model this is the case for 
the just mentioned schemes. The preferred choice may be motivated by the physical problem. 
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Since in the so defined 2PI approximation one has a functional structure solely defined in terms 
of the approximated one- and two-point functions ip and G, path integral properties and relations 
such as (16) do not need to hold true any longer. Thus irrespective of the chosen truncation scheme 
the n-point functions defined by the approximated 2PI functional generally may no longer coincide 
with the corresponding 1PI vertex-functions and thus WTIs derived from the 1PI formalism may 
be violated for those self-consistent approximations. 

The reason is simply seen in diagrammatic terms. The Schwinger-Dyson resummation creates 
a very restricted subset of 1PI diagrams resulting from the iterative insertion of self-energy pieces. 
In particular this implies that already the crossing symmetry is violated at orders of the respective 
expansion parameter beyond those included in ^apprx- Thus it has to be expected that in general for 
the self-consistent approximation scheme ( |l7|) the WTIs for the vertex functions derived from the 
1PI formalism are not fulfilled for the self-consistent approximations of the self-energy and higher 
n-point functions derived from the 2PI functional. For example the symmetries may already be 
violated for the self-energy E. Especially for systems in the Nambu-Goldstone phase this implies 
that the self-consistent propagators may not comply with Goldstone's theorem. 

In order to cure this problem we supplement the 2PI approximation scheme ([l7]) by an additional 
effective action r apprx M defined with respect to the self-consistent solution as 

A apprx Lit J — 1 apprx [tp, G[ip]\ with G[ip] defined by 

^Papprx^, G] 



5G 



< 19 » 



G=G[<p] 



Here G[ip) is given as the Schwinger-Dyson solution in presence of a given mean field <p. Strategies 
like this date back to Baym and Kadanoff JTH|, the equivalence of both functionals T[cp, G] and 
T[<p] at the exact level was shown in ref. ||, while the consequences for truncation schemes were 
discussed, e.g., in jl(]] in the context of a background field formulation. Here we will show that 
independent of the chosen trancation scheme r apprx [v?] permits to construct proper vertex functions 
on top of the self-consistent solutions of fll7[), which then obey the symmetries as described by 
the usual 1PI WTIs. In diagrammatic terms, to any ^-derivable approximation this procedure 
supplements that minimal set of diagrams which is needed to recover the symmetries. Effectively 
can be considered as a background field and r apprx [c/j] permits to construct the linear response of 
the system with respect to a fluctuation in ip around the self-consistent solution of the equations 
of motion (17). The latter coincides with the stationary point dp of the effective action functional 

( ^r apprx [^, G] f 5T a , pprx [<p,G}5G{2['p]\ 



r 

*- apprx 



apprx [ 



Sip 



tp=tp 



12 ' G=G[<p],<p=<p (20) 



<5r apP rx [tp, G] 
dip 



G=G[tp],tp=tp 



i.e., ip and G[<p] solve ([H 

Now since the derivation of the WTIs (15) only relies on the symmetry of the r apprx [<£>, G\- 
functional under the 0(A r )-transformations of <p and G as a contravariant vector and a contravariant 
tensor of 2 nd rank, respectively, and this by construction holds true for the approximated functional, 
relation ( |i~5|) is also valid for the approximation. Therefore we immediately conclude that the non- 
perturbative effective action functional ( |l9| ) is O(N) symmetric in the sense of the classical action, 
i.e., 

d(l) 5Ta ^; m (ry f p{=0. (21) 
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Thus, if we define self-energies and higher proper vertex functions by the usual definition as multiple 
derivatives of T, then these functions fulfill the usual 1PI WTIs, and thus obey all symmetry 
properties derived from them. We call the so constructed functions external, as they do not take 
part in the self-consistent scheme but rather are calculated as a function of the frozen self-consistent 
solutions of (|l7|), the latter being referred to as internal. While external and internal quantities 
coincide at the exact level, this is commonly no longer the case for the here discussed approximation 
schemes based on a truncated ^-functional. Below we drop the label "apprx" . 



A. Goldstone's Theorem 



As a first quantity we discuss the external propagator defined from T by 



(^ext 



ext)lj,2k 



5(f[5cp% 



<p=<p 



which fulfills the WTI 



(22) 



(23) 



as can be immediately seen by taking the functional derivative of (|2l|) with respect to tp\ an d using 
the fact that T is stationary for <p = (p. For a translationally invariant situation (for instance for 
the vacuum or thermal equilibrium) this can easily be Fourier transformed leading to 



(Oi*(P = 0)(r a )W 



-(M 2 ) jfc (T%V' = 0, 



(24) 



where M 2 denotes the (thermal) mass matrix for the field degrees of freedom. This reflects Gold- 
stone's theorem: There are as many massless states as group generators which do not annihilate 
the mean field, i.e., the number of massless Goldstone bosons is the dimension of the symmetry 
group minus the dimension of the symmetry group of the mean field. In our case the symmetry 
group is O(N) and the symmetry group of the mean field is 0(iV — 1), which means that we have 
(N — l)N/2 — (N — 2)(N — l)/2 = N — 1 massless Goldstone Bosons in the spontaneously broken 
symmetry phase. 



B. The external propagator 

Now we explicitly construct the external propagator. From ([19]) we find 

. 5 2 T[^G] 6G{%' 



lj,2k 



S 2 T[ip,G] 



i^ + /d(3'40 3 
{5^2 Jc 5ip[i 



(25) 



The derivative of G with respect to the background field (p can be expressed through the identity 



jfd(2') (G -1 )ij,2'fc'^2 = Siffij- 



Taking its derivative with respect to the mean field yields 



d(2') 



— (j "2'2 + \y )lj,2'k'—^—[ 



(26) 



(27) 
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implying 



5G\ 



12 



/ d(l'2') S 3j^&&&*. 

Jc StpL 



S(p l 3 jc 

The functional derivative of G~ l , which is a three-point function 



(28) 



A 



(3) 



^j,2k;3l- ^ . 

can be implicitly expressed in form of a Bethe-Salpeter equation (BS) 



A 0) _ r (3) _ ■ 

2 Hj,2k;3l ~ 1 lj,2k;3l 1 



/ A(^A'^."A"\ r^ 4 ) fil'l" frm'm" k($) 

J c ^ ' lj,2fc;3'Z',4'm' ( - T 3'3" ( - T 4'4" A 3"Z",4' 



m";3r 



(29) 



(30) 



We denote the resummed, i.e., two-particle reducible n-point vertex functions by A( n ) while the 
corresponding irreducible parts are denoted by . Constructed through <3? the latter 



r 



(3) 

lj,2k;3l 



r (4) _ _ 2 

1 lj,2k;3l,4m ~ 



G=G[^],¥)=<^ 



(31) 



G=G[tp],<p=<p 



are functionals of the self-consistent propagator G. 



In terms of diagrams the BS-equation (30) can be depicted as follows: 



iA< 3 ) 





(32) 



The different external points of the vertex functions, dots and stars, indicate the different kind of 



points, which are separated by a semicolon in (30). This specialty of those diagrams clearly shows 
that only certain channels are resummed namely those which are not contained in the self-consistent 
scheme. As explained in [I] the Dyson resummation implies s-channel ladder resummations, while 
the here given BS-equation creates t- and n-channel ladders, this way recovering crossing and all 
global symmetries. This again shows the virtue of the ^-functional method. Through the effective 
action functional T[<p] it permits the construction of fully symmetry preserving vertex functions void 
of double counting, although each dynamical equation, the Dyson and BS-equation, for themselves 
violate these symmetries. The kernel of the BS-equation is symmetric with respect to the 
simultaneous change of the pairs 1,2 with 3,4 and 2PI for separating these pairs. In addition the 
^-functional ensures that especially the crossing symmetry is also recovered for the counter terms, 
a property necessary to renormalize the proper n-point functions. 2 This is non-trivial since the 
counter terms appear to arbitrary order in the sense of an h- or coupling constant expansion. 



2 Power counting and Weinberg's convergence theorem show that in our case the superficial degree of divergence of 
a diagram 7 is 5(7) = 4 — E, where E is the number of external legs of the diagram. This means that due to the 
BPHZ-renormalization procedure only proper n-point functions with n < 4 have to be renormalized. 
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Using (|30|), (|3~l|), ( |A15 ), and (A4) in (25) we finally obtain for the external self-energy 



(S ex t)lj,2fc 



- / d(l'2') - 

2 Jc (J^^I^V'i'^^' 



f' c fc ,^l'2' + c J, k 



<W¥>2 

1 /" aMaW'a!'\ r (3) r fj" r k'k" A3) 

~2 / ^ ' 3'j',4'fe';lj (jr 3'3" (jr 4'4" ii -3"j" 

Again this can be easily expressed with help of diagrams 

Q 



4"fc";2fc 



(33) 



G=G[ip],ip=ip 





(34) 



where the last term of the external self-energy, the BS-part £bs > needs special care with respect to 
renormalization. The external self-energy S ext , here derived for a general ^-derivable truncation 
scheme is frequently called the mass matrix, as it is used to determine the masses of the fluctuations 
in the phase of broken symmetry. 



C. Solution and renormalization of the Bethe-Salpeter equation 



The BS-equation now given in short hand four-point function notation 

-1G13G24 



A (3) = r (3) + r (4) G (2) A (3) ) where ^2)^ 



can be solved in terms of the four-point function A^ 4 ) as 

A (3) = r (3) + A (4) G (2) r (3) ) 

where the four-point function A^ 4 ^ is defined through two equivalent BS-equations 

a( 4 ) = r( 4 )+r( 4 )G( 2 )A( 4 ) 
= r( 4 )+A( 4 )G< 2 )r( 4 ). 



(35) 



(36) 



(37) 
(38) 



The latter four-point BS-equation is identical to that already considered in the first paper [I] of 
this series, however with two differences: In I only s-channel resummations were considered which 
restricts the momentum arguments to forward scattering, and secondly the BS-equation was used 
solely at the vacuum level. Here, however, the BS-equation truly acts at finite temperature. The 
latter implies that for the renormalization procedure vacuum and finite temperature pieces need to 
be separated along the lines given in [I]. Again however the only subdiagrams to be renormalized 
are the four-point functions which lead to the renormalized vacuum function A defined in [I]. 
Considering the difference between the four-point BS-equation at finite T and that at vacuum 
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(39) 



leads to the following expression 3 

a( 4 ) _ A( 4 ; vac ) = r (4) - r (4;vac) + r (4) G (2) A (4) - r( 4; ^)G( 2,vac )A( 4;vac ) 
= r( 4 ) - r( 4;vac ) + A (4) G (2) (r (4) - r (4;vac) ) 

, ^(4)gi(2)p(4;vac) _ p(4;vac)g<(2,vac) ^(4;vac) 

= p( 4 ) - r (4;vac) + A (4) G (2) (r (4) - r (4;vac) ) 

_|_ (r( 4 ) - r(4;vac)-)£<(2,vac)^(4;vac) 

+ A^G^fA^ 4 '™^ - r( 4 ; vac )c( 2 > vac )A.( 4;vac )) 

V v ' 

p(4;vac) 

- (A (4) - A( 4 )G( 2 )r (4) )G (2 ' vac) A (4;vac) 
r< 4 ) 

= p( 4 ) - r (4;vac) + A (4) G (2) (r (4) - r (4;vac) ) 

_|_ (r( 4 ) - r(4;vac)-jg<(2,vac)^(4;vac) 

+ A (4) (G (2) - c( 2 . vac ))A( 4;vac ) 

+ A^G^fT^ - r( 4 ; vac ))c( 2 . vac )A( 4;vac ) 

Since the I^ 4 ) -functions are 2PI when cutting the corresponding diagrams such that the two external 
point pairs become separated, the subtractions cause the explicite loops in the final expression to 
be finite such that only the vacuum quantities are to be renormalized. 

This vacuum renormalization is obtained with the same techniques as explained in [I]. In the 
case of a spontaneously broken symmetry only the renormalization description has to be changed, 
since the on-shell scheme described in [I] would lead to artificial infrared singularities. Rather here a 
"mass-independent renormalization scheme" as the MS or MS scheme in dimensionally regularized 
perturbation theory is required. As generally shown in [I], for self-consistent schemes BPHZ-like 
renormalizations are more convenient to use, especially in numerical simulations as presented in 
[II]. Thus we introduce a mass renormalization scale fl, i.e., for the counter terms we set the mass 
parameter m 2 to fl 2 . Since the UV divergences are ruled by the asymptotic behavior of the Green's 
functions at large loop momenta, the same counter terms also render the integrals finite for the 
spontaneously broken phase where m 2 = —rh 2 < 0. 

Then we can apply the same renormalization conditions as in [I] for the symmetric phase now 
however taken at m 2 = fl 2 : 

E( vac )( p 2 = 0;m 2 = /i 2 ) = 0, 
d p2 £( vac V;m 2 = /2% 2=0 = 0, (40) 
A (4;vac) (s = £ = u = 0; m 2 = fj 2 ) = T^™ c \s = t = u = 0; m 2 = fi 2 ) = ^. 

Herein s, t and u are the usual Mandelstam variables for two-particle scattering kinematics. This 
procedure shows that those counter-terms necessary to render the self-consistent Dyson-equation of 
motion finite also enter here for the BS-equation. Especially the crossing symmetry of counterterms 
is also recovered consistently with the chosen approximation for <]?. 

Given the renormalized A( 4;vac ), Eq. p^ ) constitutes a regular integral equation to determine 
the finite T-dependent function A^ 4 ). The steps towards A^ 3 ) and finally S BS only involve to close 



3 As explained in detail in [I], in this subtraction technique all vacuum functions are contour diagonal, i.e., they 
vanish for arguments with mixed vertex placement 
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a further loop which is rendered UV-finite considering the respective differences 

a( 3 ) — A( 3;vac ) =r® — r^ 3;vac ) + A^G^r^ 3 ^ — A( 4;vac )G( 2 ' vac - ) r( 3;vac ' ) 

V v ' 

UV-finite (Al\ 

s BS - =r^ 3 ^G , ^ 2 ^A < - 3 ^ — r^ 3;vac ^G ( - 2;vac ^A( 3;vac ) . 

UV-finite 

Here it is important to notice that in purely scalar theories without derivative couplings three- 
point vertices are only logarithmically divergent. Further due to field-reflection symmetry of a 
theory with no generic three-point interaction term in the classical action no such counter term is 
necessary to render the three-point functions I^ 3 ) and A^ 3 ) finite; only counter terms for four-point 
functions appear in the renormalization scheme cf. (|4l]). 

This completes the proof that also for the non-perturbative higher order vertex functions re- 
sulting from r the renormalization parts can entirely be defined at the vacuum level. 



D. Comments 

The main result of our symmetry analysis can be summarized as follows: Truncated expansions 
of the 2PI functional T[<p, G] with respect to parameters, like h or 1/N, consistent with the lin- 
early realized symmetries yield approximate functionals which are symmetric, when mean fields 
and Green's functions are transformed as 0(A^)-vectors and 2 nd -rank 0(A r )-tensors, respectively. 
However, the solutions of the dynamical equations of motion may imply a resummation to arbitrary 
orders in the considered expansion parameter 4 . Then this resummation is incomplete and crossing 
and 0(A^)-symmetry may be violated at orders beyond those explicitely included in the chosen 
approximant of <&. This includes a violation of the (3- function of the running coupling constant at 
the very same order (see also [17]] and [II]). 

The techniques developed in [I] can be used to renormalize the self-consistent self-energy and 
the 2PI T-functional with temperature-independent counter terms. As usual in the here considered 
case of a spontaneously broken symmetry the T-functional can be rendered finite at a mass scale 
I? > 0, i.e., in the symmetric Wigner-Weyl phase. Thus no particular problems with massless 
degrees of freedom appearing in the Nambu-Goldstone mode arise, while symmetries are violated 
at the self-consistent level. 

The symmetry is appropriately repaired through the non-perturbative approximation ( |i~9| ) for 
the 1PI effective action functional. It was shown that all symmetries are recovered in the sense 
that the approximate proper vertex functions defined from ( |20[ ) in the usual sense, i.e., by 

f, „ = i/^L m 

dipi--- dipn 

fulfill the usual WTIs of the proper vertex-functions and can be renormalized with temperature 



independent counter terms by the techniques developed in [I] and Sect. QIC. 

Based on the symmetry violating intermediate Dyson resummation, however, it is expected 
that certain relics of that internal stage are still present also in these symmetric proper vertex 
functions 5 . Since the propagators, which define the kernels of Bethe-Salpeter equations, are given 



Exceptions are tr unca tions right at zero order as all zero-order terms are self-replicative, see the 1/iV-expansion 



VC 



discussed in sect. 

5 This is essentially true for all approximations other then self-consistent ones, e.g., in perturbation theory the 
internal structure is given by free particle properties! 
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by the symmetry violating self-consistent scheme, already the threshold structure is expected to 
deviate from the correct behavior, since the Nambu-Goldstone modes falsely appear with finite 
masses for these internal lines. 

One could think to include the external propagator into the self-consistent scheme, i.e., to 
consider all quantities as a function of G ex t rather than the Dyson G. This could indeed further 
improve the approximation. However, since implicitly this defines a new <]?-derivable scheme, with 
a $ defined by the 2PI diagrams leading to S ext , it has the consequence that the symmetries may 
again be violated, although at a much more "remote" level. 



IV. APPROXIMATIONS OF THE ^-FUNCTIONAL 

In this section we first give a functional derivation of the ^-functional with the techniques 
developed in Sect. ^ up to order h 2 to exemplify the reduction to 2PI diagrams and the symmetry 
properties of the approximations for the functional. As an application we give a diagrammatic 
derivation for the functional up to order A 2 . 

We write the Lagrangian of the 0(iV)-model in the following form: 

^ = \{d^){d^) + ^) 2 - (43) 

It is immediately clear that in this case for both, the loop expansion and the coupling constant 
expansion, the corresponding approximations for the ^-functional are symmetric under 0(N)- 
transformations, where ip and G transform as a vector and a tensor of 2 nd rank, respectively. 



A. The ^-functional up to order ti 2 

According to our discussion in Appendix [A] we have to calculate the functional 



Zi[J'] =Aexp ( -S[(f/ , J'] ) I D$"exp 



iA 



d(i2) ($r%-,2*(^W')2 



- ^s jklm J d(i) [K^'tifW'WT + ^^''){{fW')W n i 

where Sjkim = $jk&lm + bjibkm + Sj m 8ki 



(44) 



and where we have explicitly reintroduced the internal field indices. The ^-expansion can now be 
generated following the standard technique for the usual perturbative calculation (see, e.g., fl8| ). 
All path integrals of the h expansion can be found be taking functional derivatives of the Gaussian 
path integral 



Z 10 [J',K] =N J Dd>"exp (7d(12) {^i^'W't 

ik 



exp 



iTrln^-VM 2 )-^ 



d(i2) ygidjKu 



(45) 



with respect to the new independent local auxiliary source K which has to be set to at the end 
of the calculation. We have to substitute only 
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in the polynomial expression for the expansion of exp(Sj[Vh(f)" , <p]). 
In our case this leads to 

W X [J'\ = -ihln[Z 10 [J' , K = 0](1 + Hz? + h 2 zf ] )} + 0{h z ) 

= S[^ ,J'] + ^TrH^ l M 2 ) (47) 

+ h 2 (^6 jMm J d(l) + y<Wi'^'m' jT d(12) A^n m '^2 l '4}] ■ 

Now we need to substitute <p' = instead of 9? in this expression. To that end we have to find the 
expression for 

up to order h. Since ipo is the stationary point of S[(f>', J'] at fixed J' we have 

Wo, J' 



Further from the equation of motion for (p' we have 



SPo- ( 4 9) 



2A- 



<5J' 



_ =*12 ' 
12 



From this we get 

ft 5 



2 <5(J 

and thus 



^Trln^/M 2 ) =-Trln(G- 1 M 2 ) 



\\ 2 5 jklm 5 jim , m , J d(12) ^G^Gr 2 m 'G^V 2 ' + 0(/i 3 



and 

^2 



(50) 



— Trln^/M 2 ) = ^X5 klmn [ d(l') Gft^Gft + O(tf) (51) 
')i 2 



(52) 



% ] = V m <V fcW jf d(12) ^i^G^'G^VC + 0(h 3 ) (53) 

Gathering all terms finally yields 

T[<p, G] = - I / d(12) ®£{G 12 - ® 12 ) + ^[J ] 

= S[<p] + - Tr HG^/M 2 ) + — «5 jWrn jf d(l) Gjfci? (54) 

+ / d(12)^ 1 Gt 2 fc 'G^Gr 2 m VC + 0(/ i 3 ). 

J c 

This is indeed invariant under global O(iV) operations if the mean fields tp and G are transformed 
as tensors of 1 st and 2 nd rank respectively. 
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B. The ^-functional to order A 2 

Now we like to give an example for the application of the diagram rules for the ^-functional 
and the equations of motion for the self-consistent field and self-energy. 

Prom flAlOD we see that the ^-functional is given by diagrams which obey the Feynman rules 



with the point vertices of a field theory with an interaction Lagrangian given by Si [(f), <p]. The Lines, 
connecting these point vertices denote full propagators \G rather than perturbative propagators. 
As argued in the paragraph after equation (|A17 ) one has to keep only two-particle irreducible 



closed diagrams with at least two loops for i<£. 

For the 0(iV)-model this leads to the following diagram expression for i$ up to order A 2 in the 
coupling-constant expansion scheme 



i<%,G] = X X ( ) X <TZ>- (55) 




The corresponding analytic expression reads 

h 2 x 



<%,G] =—8 ]klm d(l)G^T 
J c 

+ J d(12) <AGuC&GT'4 (56) 

+ l ^8 jklm 8 rm , ml J d(12) G$G$G&G& m ' + 0(A 3 ). 
The equations of motion for the approximation are given by ( A16- Al~7| ). For the mean field we find 



+ — g — 8jklm.Sj'k'l'm' J d(2) G12 G l i 2 G™™ tp J 2 — 0. 
In terms of diagrams this equation of motion reads 



(57) 



-i(-D + m>i = + 

The equation for the self-consistent self-energy follows from (Q) 

— i£lj',2fc = —^-fijklrnG^S^ + ^ 8jpqr$kp'q'r'¥lGl2 ^12 ^2 

^ x a r>PP' nil' n rT> 

q — jpqrO hp' q'r ,(j 12 Lx 12 Lt 12 (59) 

According to (^1-34) the diagrammatical elements for the calculation of the external self-energy 
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are given by 



di$>[<p,G] 



ir (3) = X. C ) 




(60) 



srw= V f) (62) 



V. A SIMPLE EXAMPLE: THE HARTREE (PLUS EXCHANGE) APPROXIMATION 

To illustrate the above formal considerations in this section we numerically solve the equations of 
motion for both the self-consistent and the symmetric effective self-energies for the approximation 
keeping only terms to linear order in the explicitely appearing coupling A in all above equations. 
This defines the Hartree approximation including the corresponding bosonic exchange terms 6 . 
Especially c.f. (^), becomes 

*[<P,G\ = ±6 jklm [ d(l) G&Gft, (63) 
J c 

where we have set h = 1. 



A. The vacuum case 

In the following we can restrict ourselves to time ordered functions, so that in this section all 
propagators stand for { }-propagators and all vertices for { — }-vertices. 

First we have to find and solve the self-consistent equations of motion (||) for the approximation 
given by the $ functional (|6^). For a given mean field (p\ we find: 



\j,2k 



-n 1 + rh 2 )5 jk - ^(2<pij(pi k + (p\5jk) 



5 {d) (x 1 -x 2 ). (64) 



Since the vacuum is homogeneous in space and time the mean field is constant and pointing in one 
arbitrarily chosen direction. Therefore it is convenient to express all vector and tensor quantities 
in terms of their components parallel and perpendicular to (p. The perpendicular modes are those 
of the Goldstone bosons (e.g., the pions in the linear sigma model). The corresponding projectors 
are 

pjk _ §jk _ 
pjk _ 



6 The latter are of next-to-leading order in the semi-classical 1/N expansion. 
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In momentum representation the self-energy components derived from the approximation to <3? 
given by (|63[) become 



d d l 



Correspondingly the Dyson equations (|17j ) decouple with 

G±(p) 



[(N + l)G ± (l) + G«(l)], 



(27r) dl 

^2[(^-l)G±(i)+3G|| 



(66) 



Gii(p) 



" p 2 - M| + ir/ ' 
1 

p 2 - M 2 + i?7 



with M]_ = - m 2 + S_l, 
Mf = y^-m 2 + S||. 
With these definitions the equation of motion for the mean field reads 

r2 



(67) 



<pi(M$ - Xif ) = 0, 



(68) 



which shows that either the mean field or the bracket vanishes. In the first case, the symmetric 
Wigner-Weyl mode, in the second the Nambu-Goldstone mode is realized. In our case of negative 
—m 2 the spontaneously broken phase is realized in the vacuum. 

In the following we shall use the convention and renormalization prescription according to (PT 



In this renormalization scheme the self-consistent gap equations (66-67) become 



Ml = e ± (M ± ,Mu 



^ ->2 ~ 2 , ^ 

—op — m H 7? 

2 * 32vr 2 



(N + 2) ft 2 - (N + l)Aff - Mf 



M 2 

+ Ml(N + 1) In ( ) + M 2 In 



M 2 = e||(AfjL, Af||) := 



2 Y 32vr 2 



(JV + 2)/2 2 - (N - l)M'i - 3M, 



r2 



+ Af£(iV- l)ln 



^2" 



Mi, 



(69) 



The equations show clearly that the physical results are independent of the choice of the mass 
renormalization scale ft since a different choice of ft can be compensated by a (finite) renormaliza- 
tion of the coupling A and the mass parameter in. For more general ^-derivable approximations, 
which contain real two-point contributions to the self-energy, also a wave function renormalization 
is needed (see @ ) . The introduction of the mass renormalization scale ft is necessary only because 
we need a "mass independent renormalization scheme" [IS, 2C] which uses the fact that a the- 
ory with symmetries can be renormalized with counter terms of its symmetric realization in the 
Wigner-Weyl mode also for the spontaneously broken Nambu-Goldstone phase. In addition mass 
independent renormalization schemes have the advantage to avoid "renormalization induced" in- 
frared divergences which would appear with the massless modes within an on-shell renormalization 



18 



scheme. It is also clear that for special choices of jl our class of renormalization schemes contains 
the minimal subtraction scheme (MS) and the modified minimal subtraction scheme (MS). 

Prom ( |69|) it becomes immediately clear that by the self-consistent ^-derivable approximation 
the symmetry is explicitly broken in the Nambu-Goldstone mode of the theory, i.e., for m 2 = 
—rh 2 < which provides ip 2 > 0: The transverse modes should be massless according to Goldstone's 
theorem which does not hold true according to (|69|). 



As we have seen in section III we can define an approximation to the self-energy respecting 



the underlying symmetry by (25). For that purpose we need to solve the equation of motion 

e kernels, here given in the momentum represi 

^jk;l,m = ~~ -^{SjkSlm + ^jl^km + 8j m 5kl) ■= — — Sjklrrn 



From fl31|) we obtain the kernels, here given in the momentum representation 

(70) 



(3) 

Fjk;l = ~H<Pk8jl + Vj&kl + tpiSjk) = —^jklmfm- 

Since due to translation invariance <p is independent of the four-momentum both and are 
constant and effectively only one-point functions. From this it is clear that A® is effectively only a 
two-point function and thus depends only on the momentum attached to its third argument. This 
becomes also immediately clear from the diagrams ( |32"| ) and (j34|). In our special case the kernels 
r( 3 ) and r( 4 ) are point vertices 

ir( 3 ) = V and irw = X (71) 



and thus from the iterative solution of the equation of motion (^TJ) and ( p3| ) for A^ 3 ) we conclude 
that the external self-energy associated with the Hartree approximation is given by the RPA bubble 
resummation 



To find the analytic expressions for this diagrammatic equation we specialize ( p0[) with the kernels 
(|70|). This leads to the equation 

= -A* iWm y> ra + / ^G l ' l "(l)G m ' m "(l - k). (73) 



According to ( |40[) for renormalization we have to subtract the expression where in (73) both 
propagators are set to free propagators with the mass parameter set to the mass renormalization 
scale jl and k = 0. 

From the symmetry properties of A^ 3 ) it follows that it must be of the form 

A$,(*0 =j (Al 3) (fc)Pf + t$\k)pf 

+ Af(k)Upf + v k P j l) (74) 



+ Af ik) (y'pf + /p,f 

Indeed after some algebraic manipulations with the propagators (|67|) we obtain the algebraic linear 
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equations of motion for the four independent components of A 3 as 

Af\k) = -A + ^ [(N + l)L ±t± (k)Af(k) + (Af (fc) + 2Af (fc))L| U |(A ; ) 

A£ 3) (fc) = -A + ^ [(AT - 1)L ±)± (A;)AS 3) (A;) + (3Af (fc) + 2Ai 3) (fc))%||(fc) 
A^(fc) = -A + AL|| i± (A:)A^ 3) (yfc), 
Ai 4) (fc) = -A + Ai| M |(A:)Af ) (fc). 
Here we used the abbreviation 

L a ,p{k) = L Pta {k) = i/ 7^4 {G a (l)G p (l — k) — [D(l)] 2 } 



(75) 



with a,/3 € {_L, ||} and L>(Z) 



(76) 



Z 2 — /2 2 + i?? 

According to (|3^) and (^) the external self-energy is given by 

(Sext)ifc(p) = s ifc (p) - ^(A® (p) - rg) 



S . fe _ ^ {[3A + A (3 )(p) + 2Af (p)]^| jfc + [A + Af(p)]P I 



(77) 

-JLjfc}- 

In order to prove Goldstone's theorem we have to set p = and use the equations of motion (68) 
and (|69|). The crucial property is 



Mi = Sj. - Sii = -A[9i (M,f) - 6i(Aff) - (Mjf - M|)e 2 (A )], (78) 



where the functions Q\ and 2 are given in appendix cf . Q and Q. Using Q together 
with 

^cxt± = ^ 2 -^ 2 + 5]cx t ±(0) (79) 



and (77) we find indeed M 2 xtl _ = as it should be according to Goldstone's theorem: The (N — 1) 



degrees of freedom perpendicular to the mean field are the Goldstone modes corresponding to the 
spontaneous breaking of the O(N) to the 0(A^ — 1) symmetry group of the mean field. 

Finally we remark that the Hartree self-energy for the case of the Wigner-Weyl phase of the 
theory, where <p = 0, is identical with the external self-energy and thus fulfills also the WTI for 
the self-energy, which indeed in this case is trivially fulfilled. 

In Fig. [1] we plot an example for the above considerations: We roughly fit the parameters of the 
self-consistent self-energy such that the properties for pions are satisfied, namely M± = 140MeV 
and My = 600MeV, cp = f n = 93MeV, and n = 4. The equations of motion (|69]) were solved for 



m 



and jl 2 . The renormalization conditions (f40|) were used. The plot of the pion mass clearly 
shows that Goldstone's theorem is recovered, since (Af[ ext) ) 2 (s = 0) = although the internal 
Green's functions do not fulfill the Ward identities. This violation of the symmetry properties by 
the self-consistent self-energies is clearly seen in the spectral function for the "cr-meson" , since its 
threshold is at yfs = = 2M± = 280MeV and not at y/s = as it should be for massless pions in 
the chiral limit. 

Like the here constructed external self-energy where due to the finite mean field (p the RPA- 
bubbles contributed to the self-energy expectation value ip(x)<p(y) ((j) 2 (x)(j) 2 (y)} also other symme- 
try preserving two-point functions are given by the very same RPA terms, and thus result from 
the linear response of the system due to fluctuations around the Hartree solution. A prominent 
example is the correlator (j fi (x)j l/ (y)} of the Noether current (0), which then is conserved. 




B. The finite temperature case 



Since the counter terms at finite temperature are the same as for the vacuum case we can 
immediately write down the renormalized equations of motion where we make use of the vacuum 
functions defined in (^) and the explicitly T-dependent finite part of the tadpole diagram: 

;\ r Ail :\ poo 

(T) (M) = | J __ 27rnT{lo)6{l 2 _ M 2 ) = _ J dwv / CJ 2_ M 2 nTM) (80) 

with the Bose-Einstein distribution function nx- 

From the fact that the self-consistent { }-propagator reads 

G—(p, M) = 1 . - 27Tm( Po )5 ri (p 2 - M 2 ), 

p z — M z + 17] 

1 ( \ 1 \ 

v (x) = Im 

Z7T1 \x — ir] x + vqj 

and that for our tadpole integral we are allowed take the limit rj — > +0 in the explicitly temperature 
dependent part we find the renormalized gap equation at finite temperature 

= p(Mj - Xip 2 ), 

Ml = 0j.(AfL, M||) + (N + 1)6 (T) (M ± ) + 6 (T) (M||), (82) 
Mf = 0||(Mj_,M||) + (N- 1)6 (t) (M_l) + 36 (T) (M||), 

where 0j_ and 0n are defined by (p9|) 7 . The solutions of these gap-equations are shown in Fig. ^ 
which shows clearly a first order phase transition behavior. There exist two "critical temperatures" : 
For < T < T c \ there exists only one solution with <p ^ 0, for T c \ < T < T C 2 two solutions with 
(p ^ and the symmetric solution with ip = while for T > T C 2 the only solution is symmetric. Also 



(81) 



Note that although the self-consistent masses M± and M« are temperature dependent the counter terms used 
to render 9i and 0|i finite are temperature-independent since the counter terms oc Ex and E|| are due to the 
subtraction of hidden vacuum-divergences of the four-point function as explained in detail in [I] . 
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Self-consistent o-mass at finite temperature 





0.7 r 




n r - 




u.o 






> 




CD 




(3 


0.3 - 





0.2 






0.1 - 









-0.1 - 



broken phase 1 
broken phase 2 
symmetric phase 



> 

CD 

(3 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
T[GeV] 



Self-consistent jr-mass at finite temperature 
0.6 



0.5 
0.4 
0.3 
0.2 
0.1 


-0.1 




' broken phase 1 
broken phase 2 
symmetric phase 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
T[GeV] 



FIG. 2: The solutions for the gap equations (p2[) 



the equations of motion for the external self-energy (|77|) remain formally the same but has to be read 
within the real-time {— +}-matrix formalism. As cited in [I] appendix A the analytic properties of 
two-point functions lead to the conclusion that the retarded (and also the advanced) propagators 
and self-energies decouple from the other degrees of freedom, so that the simple algebraic properties 
are valid for them as in the vacuum. We make use of the property 



F R (p) = ReF~(p) + itanh 



ImF— (p), 



(83) 



which holds true for any amputated two-point function. 

In our special case we can take advantage of these analytic properties also for the functions 
, since effectively these are two-point functions as well. Thus we can use the vacuum equations 
(|7q) for the retarded functions without changes. Furthermore it is clear that these functions do 
not contain any renormalization parts except the already removed pure vacuum divergences. 

Since the effective mass at p = is identical with the second derivative of the effective potential, 
defined by 



V eS [<p]6V>( P ) = -f[<p]\ 



i/3=const. 



(84) 



its value should be > for a stable solution, i.e., for a minimum of the effective potential rather 
than a maximum which provides an unstable "tachyonic" solution. The explicit calculation shows 
that the solution, denoted by "broken phase 2" in Fig. ^ are unstable. This shows that we find a 
phase transition of first order, i.e., a discontinuity in the order parameter ip 2 = M|/A (see Fig. 1). 

The effective masses for the stable spontaneously broken phase together with the spectral func- 
tion for the ex-meson are depicted in Fig. ||. 

The calculation clearly shows the symmetry violations of the underlying self-consistent Dyson 
approximation: Although the WTIs and Goldstone's theorem are fulfilled for the external propaga- 
tor remnants of their violation by the internal propagators are present: The low-energy threshold 
behavior of the "cr-meson" is not correct since the "pion" mass within the internal propagator does 
not vanish. Also the phase transition comes out to be of first instead of second order as it should 
be. 
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External o-mass at T=150 MeV (stable solution) External o-mass at T=150 MeV (stable solution) 




External n-mass at T=1 50 MeV (stable solution) External o-mass at T=1 50 MeV (stable solution) 




a-spectral function at T=150 MeV (stable solution) 



-Im G Rex, [GeV" 2 ] 




FIG. 3: The effective external masses at a temperature of 150MeV. The effective external 7r-mass indeed 
vanishes at po = p = as predicted from Goldstone's theorem. The spectral function of the a-meson shows 
that at high temperatures its strength becomes more peaked and the maximum shifted to lower momenta 
than at T = 0. 



C. Leading order large-N 

In the context of symmetries the large TV expansion scheme is of particular importance. Here N 
denotes the number of fields, e.g., in SU(./V) or O(N) theories. The counting depends on the type 
of theory and is defined such that the classical action scales like N, c.f. fl4|, gjj . As unrestricted 
loops scale like N this implies that coupling and mean fields scale like 

unrestricted loops oc N, A oc — , <p oc V~N (85) 

for the here considered O(N) model in the Nambu-Goldstone phase. As the WTIs concern the self- 
energies in the first place we shall first dicuss the 1/N expansion at this level before we comment 
on it in the context of the 2PI-functional formalism. 
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l/7V-expansion of the self-energy 

The 1/N counting scheme has the remarkable feature that the leading order (LO) leads to zero 
order terms, i.e., oc (1/iV) , in the self-energies. This implies that iterative self-energy insertions 
of zero order contribute in LO which amends the entire LO self-energy self-consistently to be 
constructed within a corresponding Dyson resummation. Thus the resulting LO-propagator is 
created fully self-consistently and at the same time accounts for all terms at LO. The latter fact 
guarantees that for the LO self-energies the symmetries, i.e., the corresponding WTIs, are fulfilled. 



Since the above derived Hartree + RPA (external) self-energy, given by Eqs. (gg) to (77) 
indeed includes all zero order terms in 1/iV-expansion, the leading order can simply by obtained 
by retaining the according large N limit terms 



_ x LO 

1s lo = J _y_ + JsQzt + AqqzIl 



+ ■■■} . (86) 



Naturally this leads to counting factors different from those in the RPA result ( |72[) as different 
contractions of the (y?</?) 2 interaction term lead to different 1 /iV-orders (see, e.g. , |l4j ) . The counting 
further implies that (a) solely loops in the Nambu-Goldstone modes survive and (b) therefore the 
bubble sum only contributes to the a-meson self-energy Sy. The explicit result is 

vL0 _wv t jfLafd), with (Mi°)V = o, 

1 ' ' (87) 



2 

:£V + 



2\2 /V rLO 
LO _ v LO ■ V A iVjL J ■ 



1 - AiVLLO ' 



Here L^°j_ denotes the pion bubble loop, c.f. (76). The result clearly shows that for the broken 
phase (pi the "pion" -mass vanishes for the zero-order propagator. The renormalization of 
Eq. (^7j) is done in the same way as described above for the Hartree case and of course also the 
renormalized "pion" -mass vanishes for both, the vacuum and the finite-temperature case. The 
WTI for the self-energy ( |87|) are fulfilled by construction, if all terms of a given order are included, 



and thus both the 0(A r )-Noether current (11) is conserved and the Goldstone theorem is fulfilled. 



For the zero-order result (pTl) all this can be proven by mere inspection. 



Please note that although the LO result (87) can be generated from a self-consistent Dyson 



resummation scheme, it is not ^-derivable, since the pion self energy pieces corresponding to the 
bubble series are of subleading order! Higher order approximations can even not be expressed 
in terms of a self-consistent scheme, since from the Dyson series the corresponding propagator 
contains all orders, while the terms contributing to the self energy are of limited order! 



1/iV-expansion of the 2PI generating functional 

In recent times the l/A r -expansion scheme also has been used to organize truncation schemes for 
the here considered 2PI generating functional. Thereby the counting results from the scaling rules 
( |85| ) , however excluding the internal structure of the self-consistent propagators from the counting 
in the diagrams of $ ([[]), c.f. refs. [[lj, This procedure leads to a 2PI generating function 
which is symmetric, a premise for the entire discussion in the papers. 

In leading order 1/N one recovers the diagrams for the Hartree approximation (|64| ) sect. ^ 
and the corresponding Dyson resummation result, however with the counting factor arising from 
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LO in 1/N. The internal self-energy resulting from the corresponding Dyson resummation is given 
by the above result (|8q), (|87|), however with the essential difference that the bubble contributions 
are dropped. The reason is that the corresponding 2PI diagrams, c.f. Fig. 5 in ref. [[13], are 
of subleading order. Nevertheless the pion mass is still zero, since it is not affected from this 
difference, while the sigma-meson self-energy is lacking the corresponding decay cuts and the 
WTIs at the correlator level are indeed violated. The external self-energy then exactly recovers the 
missing diagrams leading back to the result (p7|) . In fact for the 2PI-l/iV-expansion the Goldstone 
theorem can only be assured for the so called mass-matrix which is identical to the here considered 
external self-energy. A noteworthy side feature, though, is that the LO approximation gives the 
correct 2 nd -order phase transition. For details sec [fToj] . 

As the main theme of this paper we explicitely see here that a symmetric generating functional 
by itself does not guarantee that the solution of the equations of motion, i.e., the Dyson equation, 
preserves the symmetry. Rather the here presented functional scheme to construct external vertex 
functions is the minimal procedure to cure the corresponding symmetry defects. 

VI. CONCLUSION AND OUTLOOK 

In this paper we have analyzed Baym's ^-derivable approximations with respect to their sym- 
metry properties. It was shown that the self-consistent mean field and propagator solutions of 
such approximations in general do not fulfill the Ward-Takahashi identities of the full propaga- 
tor as derived from the usual 1PI generating functional, i.e., the effective action. The reason for 
that lies in the fact that, although the expansion of the 2pI-functional is done systematically in 
a symmetry-conserving ordering scheme (for instance in powers of some coupling constant of a 
symmetric term in the Lagrangian, the h- or the l/A^-expansion) the solutions of the equations 
of motion correspond to an incomplete resummation to any order of the expansion parameter. In 
general not even crossing symmetry is respected for the solutions beyond the order of the expansion 
parameter. 

Furthermore it could be shown, though, that for any such truncated Dyson resummation scheme 
it is always possible to define a non-perturbative lPI-effective action based on the self-consistent 
solution. This supplementary action indeed generates proper n-point vertex functions which fulfill 
the whole hierarchy of Ward-Takahashi identities. These external proper vertex functions are 
implicitly determined by closed vertex equations of Bethe-Salpeter or higher order type, where all 
ingredients are constructed from the self-consistent mean field and the self-consistent propagator 
resulting from the underlying ^-derivable Dyson scheme. The corresponding solutions exactly 
recover the crossing symmetry and at the same time the symmetries of the original classical action 
provided there are no intervening anomalies. The fact that the ^-functional formalism was used to 
determine the self-consistent mean fields and propagators does not only avoid double counting but 
also ensures the consistency of counter terms, needed for renormalization of the divergent integrals. 
This is valid for both, for the 2PI functional T[<p, G\ and its equations of motion as well as for 
the here discussed higher order vertex equations like the Bethe-Salpeter equation resulting from 
the effective quantum action T[<p\. In the simple Hartree approximation the symmetry repairing 
procedure leads to the well known Random Phase Approximation (RPA) . Similar features result for 
the leading order 1/N expansion of the 2PI functional. For approximation levels with genuine two- 
point self energies, i.e., for problems with damping, one arrives at Bethe-Salpeter equations of ladder 
type, as they were considered, e.g., in the context of Hard Thermal Loop (HTL) resummations 
or of photon production p^], which in semi-classical approximation leads to classical transport 
equations. 

Although the effective lPI-action formally fulfills all the symmetry properties of the original 
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action it now suffers other defects essentially resulting from the lack of self-consistency for the 
so constructed external multi-point functions. Their internal structure is inherited from the self- 
consistent (e.g., Hartree) propagators: Since the latter do not fulfill the Ward-Takahashi identities 
it is not guaranteed that they show the correct dynamical behavior. For instance in the case of 
the here considered 0(iV)-model in the spontaneously broken Nambu-Goldstone phase without 
explicit symmetry breaking the transverse ("pionic") degrees of freedom have a non-zero mass in 
the self-consistent propagator and thus all "pion" loops show a wrong dynamical behavior (e.g., 
the vacuum "cr" spectral function shows a threshold at s = 2Mj_ which is finite rather than zero!). 
Another consequence of this symmetry problem at the internal level is the wrong prediction of a 
first order phase transition. 

Both, the here considered 1PI and the 2PI effective actions T[(p] and T[ip, G] have the same values 
at the solution of the self-consistent solutions for the mean fields and propagators and therefore 
represent a non-perturbative approximation for the thermodynamical potential. This implies that 
the symmetry violation at the internal level may also leaves its traces in this thermodynamical 
potential. 

The derivation of approximation schemes that fulfill all symmetry properties of the underlying 
classical action and at the same time are fully self-consistent still remains as an open task. 

In a forthcoming publication we shall discuss the symmetry properties for the specially inter- 
esting case of local gauge symmetries. Here the background field method permits to establish 
an explicitly gauge-invariant non-perturbative effective action on the basis of the 2PI-action [16], 
which leads to vertex functions that formally fulfill the Ward-Takahashi identities. In this case, 
however, the symmetry violation by the intermediate 2PI approximations causes even more serious 
complications, namely the excitation of unphysical degrees of freedom for the gauge field propa- 
gators. The latter fact may imply artifacts in the external proper vertex functions which among 
others even may spoil thermodynamic consistency due to the false number of degrees of freedom 
in the internal lines. Here appropriate projection methods integrated into the Dyson scheme as 
recently suggested by us 23] may establish a suitable work-around for such problems. For a more 
general recent review on Schwinger Dyson equation approaches to non-abelian gauge theories, in 
particular applied to the low-energy properties of QCD, see [24]. 
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APPENDIX A: CALCULATION OF THE 2PI FUNCTIONAL 

In order to give a precise meaning for the approximation schemes considered we briefly derive 
the Feynman rules for the 2PI functional in terms of path integrals. The main line of arguments 
follows H . In order to perform a systematic loop expansion here we explicitely shall keep track of 
h factors. 

We shall restrict ourselves to theories with only scalar boson fields. The generalization to other 
field contents is straight forward using the very same functional integral techniques. Thus we 
assume the theory to be defined in terms of a classical action functional 



S[0= jfd(l) Sf(<k), 



(Al) 
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where (f> denotes a multiplet of scalar fields. Here and in the following we use the notation intro- 
duced in I: f c d(l . . . n) /(l, 2, . . . , n) denotes the cf-dimensional integral in the sense of dimensional 
regularization with the 0-component running along the extended Schwinger-Keldysh path C (run- 
ning from the initial time ij to a final time tf and back along the real time axis and then down to 
ti — \f3 parallel to the imaginary time axis and at the end taken ti — > — oo and tt — ► +00). Here 
this contour integral also implies the appropriate sums over charge-space indices. 
The generating functional with local and bilocal sources is defined by 



Z[J,B] = N / D</>exp 



'192 



(A2) 



where N is a indefinite normalization constant which will be chosen such that for the temperature 
T — > and J = B = the functional is normalized to 1. 

Shifting the integration variable of the path integral by an arbitrary field ip one obtains 



Z[J,B] =iVexp 



±s[<p\ + y d(i) j m + ^ y d(i2) b IWW2 



D<f>' exp 



1 

2h 



d(12) (G~ 1 )i20 / i0 2 + iSjtfM + l - d(l) 



(A3) 



Zi[J'] 



where we introduced the following abbreviations 

5 2 S& 



12 



+ B 



12 



S<pi8ip2 
S I [<l> , ,<p]=S[<l>' + <p]-S[<p] 



(^- 1 ) 12 + 5 12 , 

c <Vi 



^i<^2> 



(A4) 



J[ = Ji + 6 -p^+ [ d(l') iW- 



Now we want to chose J and -B such that ip and G are the exact mean field and the exact propagator 
respectively. These are defined with help of W[J, B] = — iMn Z[J, B] via 



SW[J,B] JW[J,B] 

<Pl = ~ , 1W-T12 ~ 



6J1 5B 12 
From the first condition we derive immediately that we have to chose J' = Jq such that 

6Z 1 [J'\ 



SJ> 



0. 



(A5) 



(A6) 



In this way the functional Z±[Jq] defines completely Z[J,B] and can be calculated approximately 



with well-known standard techniques (see, e.g., [|18[| ). An example will be given in Sect. IV, In 
this way any expansion with respect to the number of loops (powers of ft) or with respect to the 
coupling constant defines a ^-derivable approximation. Of course any other expansion scheme, 
known from usual perturbation theory, is feasible. E.g., recently in [pq ] and [14] investigations of 
the Large- A^-expansion description for Baym's functional were undertaken. 

For sake of completeness we give the derivation of Baym's ^-functional. To that end we only 
have to perform the ^-expansion up to first order, i.e., up to one-loop order in terms of diagrams. 
For this task it is convenient to introduce the new action 



S[<P';J'] 



1 



d(12) (G _1 )i20 / i</>2 + Si[4>', <p] + / d(l) J^i 



(A7) 
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so that the functional Z\, defined in ( |A3| ), reads 
Z X [J'\ = N f D^exp 



^S[c/>';J'} 



exp( -WxlJ') 



(A8) 



To obtain the /i-expansion we have to expand the functional integral around the solution of the 
classical field equations ip' , which is the stationary point of the classical action: 



5S[<p' ;J'} 







(A9) 



01 



and substitute y/hcf>" = 4>' — ip' : 



Z 1 [J'] = N exp <j -S[<p' , J'\ } I D0"exp 



i / d(i2) ^-^Wi+'-sAVnfM 



(A10) 



with the definitions 



(ST 1 ) 



lj,2k 



^01^02 



Si[Vn<f>", ip' \ = s[ip' + Vhcf>"; j'\ - s[Vh<p' ; j'} 



d(12) (^- 1 )i20i'^'. 



(All) 



Note that (Sf 1 )i2 and Si both depend on J' only implicitly via (p f and that 5^ [s/hip", <po\ = 

0[h k l 2 ] where denotes the monomial to order <p" k . By definition only the terms with k > 3 are 
different from 0. The most general renormalizable theory has only k = 3 and k = 4 contributions 
in Si and we shall restrict ourselves to this case. 

Now it is easy to extract the one-loop contribution (i.e., the O(^)-contribution) to the generating 
functional W' for connected diagrams: 



Wi[J'] = -iH 



lnN + j-S[(p' , J'} + In / Dc/>" exp 
S[cp' , J'} + | Tr In^n/M 2 ) + W^J 1 ]. 



(A12) 



Herein we have introduced an arbitrary mass scale M to avoid dimensionful quantities within 
the logarithm which takes account of the overall normalization of Z which is irrelevant for any 
physical quantity derived from it. Now according to ( |A6| ) we have to chose J' = Jq such that 



SWi/SJ' = 0, and we have 



(A13) 



so that we can substitute iff instead of ip' in ( A12) leading only to a modification of the functional 
W'2 to order 0(h 2 ) while the 0(7i)-part remains unchanged. 

According to ( A6) we have to chose J' = Jq such that ip' = to obtain the original W- functional. 
We also note that for this choice £f becomes G according to (|A3|) : 



W[J,B] = -ih\nZ[J,B] = S[<p] +J d(l) Jifi + \ J^d{12) B 12 ^<p 2 + W X [J' Q ]. (A14) 



28 



Now we define the 2PI effective action by the double Legendre transformation of W with respect 
to J and B. Using (|A5|) this leads to 



T[<p,G\ = W[J,B] - f d(l) Jwi -\\ d(12) (^1^ + ^12)512 

= S[p] + | Tr \n{G- l /M 2 ) + | J d(12) 12(^12 - ^12) + G]. 



(A15) 



From the derivation we note that <&[<p, G] = 0(h 2 ), i.e., in the language of diagrams it contains 
only diagrams with at least two loops. The main difference to the usual 1PI effective action is that 
the lines appearing in the diagrams symbolize full propagators iG rather than perturbative ones. 
The equations of motion are now given by the fact that we like to study the theory for vanishing 



auxiliary sources J and B. From the Legendre transformation ( A14 ) we can immediately express 
this in terms of the functional T: 



ST[<P,G\ 

Sip 

5G 



-Ji - J d(2) B 12 <p 12 = 

■2 512 = °- 



(A16) 



Using the /i-expansion ( A15) the last line reads: 



12 



G12 



2i 6$[tp,G\ 



h 5G 



'12- 



12 



(A17) 



It is clear that E12 is the exact self-energy expressed in terms of the exact connected propagator 
G and the exact mean field ip and thus ( A17|) is the full self-consistent Dyson equation. This 
implies that no propagator line in the diagrams must contain any self-energy insertion, because 
these lines denote already the full propagator. Thus G] consists of all closed diagrams with 
point vertices from the action Si[ip, 0] with at least two loops 8 , which have the additional property 
that it is impossible to disjoin them by cutting only one or two lines, i.e., all diagrams contained 
on <I> must be two-particle irreducible (2PI). This must hold true because taking the derivative 



of <I> with respect to G according to ( A17 ) defines the full proper self-energy which must be 
1PI and no propagator line should contain any self-energy insertion. Now taking this derivatives 
diagrammatically means to open any propagator line of all diagrams contained in <E> and taking 
the sum of the so obtained skeleton self-energy diagrams. This indeed implies the 2PI property for 
the diagrams defining since otherwise one could disjoin a diagram by cutting two lines, leading 
to a 1PR diagram contributing to the self-energy which by definition cannot be contained in the 
diagram expansion of the self-energy. 



In Sect. IV on the analytical example for the 0(-/V)-linear <r-model up to order h 2 we see that 
the double Legendre transformation ( |A15| ) indeed leads to 2PI diagrams for $[<p, G]. In this section 
also an example of diagrammatical derivations for the ^-functional for the same theory is given. 



APPENDIX B: SYMMETRIES OF THE CLASSICAL ACTION 



For sake of completeness we summarize the derivation of Noether's theorem [26] for symmetries 
of a classical action functional. For sake of simplicity we shall restrict ourselves again to a multiplet 
of scalar fields. The extension to more general cases is straight forward. 



Note that <j) denotes the "quantum field" , integrated over within the path integral and ip the "classical background 
field" . Si consists of the sum over all monomials with at least three quantum fields. 
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We investigate the behavior of the classical action functional S[<f>] under a general infinitesimal 
transformation of the form 



x 'f = x tl + 5x», $'{x') = $(x) + 5${x), 



(Bl) 



where 5(f) may depend on both, the fields and the space-time argument. The action functional is 
defined to be symmetric under the transformation (Bl) if its variation 



SS[4>] 



dVjzf ((f)', d'^', x') - I d d x^((f),d^(f),x) = 0. 



(B2) 



vanishes identically, i.e., without any restrictions on the fields. To derive explicit conditions for «£? 
to fulfill the symmetry condition (B2) we have to rewrite the first integral in terms of x, where we 
have to take into account the Jacobian of the volume element. In linear order of the variation we 
have 



det 



(B3) 



The transformation and differentiation with respect to the space-time arguments does not commute: 

5(d^) := d'^'(x') - = d„{s<f,) - (d^5x v )d v (f). (m) 

After some algebraic manipulations we find 

5S[4>] r . 



5S[(f)} 



d n x- 



[d v (f))5x v \ = 0. 



(B5) 



Let now 5r] a be the independent parameters of the Lie group acting on the fields and the space-time 

(B6) 



variables, i.e., (|Bl| ) reads 

5(f)(x) = r a (x,(f>)5r, a , 5x» = -fZ(x,4>)5 V a . 
This means that the symmetry condition flB5|) reads 

5S\ 



d n x 



5(f>(x) 



{r a (x, <f>) + [8 u (p{x)]f v a (x, (f))} 5r, a = 0. 



(B7) 



Since this must hold for any field configuration for which the action is well defined and the 5n a 
are independent generators of the group operation for each a there must exist a current j a , the 
Noether currents corresponding to the symmetry group, such that 



5S[(f)} 



{T a (x,(f>) + [d u (f){xW a {x,(t))} = d^. 



(B8) 



Now the classical field equations of motion are given by the stationarity of the action. Thus for 
the solutions of the equations of motion the Noether- currents are conserved. 

To find the explicit expression for the Noether-currents we go back to (|B5| ) and express it in 
terms of the Lagrangian. After some calculations we obtain 



flfjSf 5x 



'd{d^) v ) d(d^) 
This means that there must exist local functionals Sla((f),x) such that 

SSf + Sfd^Sx" = d^ a 5n a , 



+ 5^ + ^djx^. 



(B9) 



(BIO) 
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which leads to the desired explicit expression for the conserved Noether current: 

~ - **) - + W ' (BU) 

It should be noted that within special relativity only the total conserved quantities are physical 
observables, not the local currents themselves. Those follow from the continuity equation d^j^ = 
by integration over a space-like closed hypersurface of space-time. Taking the special hypersurface 
given in a Lorentz reference frame by x° = to an d x° = t\ we find 

/ j3- a -0/-JX _ / / j3: 



Q(h)-Q(t )= dx u / d 6 xd f(x)= / d J £divj = 0, (B12) 

Jt JR 3 Jt JR 3 

where we have made use of the continuity equation. Eq. ( |B12| ) tells us that the total Noether 
charge 

Q(t) = [ d 3 xj° = const. (B13) 



For sake of completeness it should be mentioned that the currents (Bll) are not determined 
uniquely by the symmetries since changing it according to 

(j'Y = d p to pfl with to pfl = -lu pp (B14) 

does not change the Noether charge ( pi3| ), and for j' the continuity equation holds true as well 
as for j. This notion is important especially in the context of particles of higher spin where the 
freedom of choice of the Noether current, in QB14 ) parameterized by the antisymmetric tensor uj p ^, 
can be used to give gauge invariant definitions of the energy momentum tensor. In general its 
canonical version 

r r)f£ 

d(d v 4>) 

is not gauge invariant. It is clear from ( Bll|) that ( |B15 ) is the "Noether current" of the symmetry 
of the physical laws against space-time translations which leads to energy-momentum conservation: 



d_ 

dt 



d A x Q Uu (x) = 0. (B16) 



APPENDIX C: SOME FEYNMAN INTEGRALS 



In this appendix we give three dimensionally regularized Feynman integrals, needed for the ap- 
plication to the tadpole approximation and the corresponding RPA-summed external self-energies. 



The techniques to obtain them can be found in standard textbooks, for instance in [27]. As usual 
we set d = 4 — 2e for the space-time dimension in the sense of dimensional regularization, /i for the 
regularization scale; 7 ~ 0.577 stands for the Euler-Mascheroni constant. 

The first integral is the vacuum tadpole with a free propagator for a mass m: 



Gi(m 2 ) 



2c 



m 



(2vr) 4 I 2 - m 2 + b] 16vr 2 



1 , / m 2 



+ 0(e). 



(CI) 
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The next integral is used to renormalize the divergent four-point sub-diagrams contained in the 
tadpole approximation as a hidden divergence: 



6 2 (m 2 ) = i 



d 4 Z 



2c 



1 



(2?r) 4 (Z 2 - m 2 + ir/) 2 2m 



d m Qi(m 2 ) 



16vr 2 



1 



+ 7 + ln 



m 



+ 0(e). 



(C2) 



Further for the calculation of the external self-energy we need the two-point function 



Lmi,m 2 (p 2 ) — i J 



mu 



2c 



16vr 2 



(27r) 4 (I 2 -m\ + ir/) [(I - p) 2 - m| + ir/] 

/ 1 \{p 2 ,m 2 ,m 2 ) 
2 + 7 + - 



artanh 



V 

9 9 9 

m\ — mj + p 
X(p 2 , m 2 ,m 2 ) 



artanh 



2 2 2 

X(p 2 ,m 2 ,m?,) 



+ 



+ 



m\ — m\ (m\\ /mim 2 \ 



m 2 



V 47T/X 2 y 



+ 0(e), 



(C3) 



where the Kallen function reads 



X(p 2 , m 2 ,m 2 ) = ^/[p 2 - (mi + m 2 ) 2 ][p 2 - (mi - m 2 ) 2 ]. 



(C4) 



For the proof of Goldstone's theorem for the external propagator we need this function at p = 
which can be expressed with help of the already defined function 0i: 



Lnii ,m,2 (0) 



1 



2 2 

m\ — m-2 



[ei(m 2 )-Gi(m 2 



(C5) 



Their expressions for equal masses read 



■Lm,m (p) 



16vr 2 



- - 2 + 7 + In 

e 



m 



A(p 2 ,m 2 ,m 2 ) 
+ 2 = artanh 



\{p 2 , m 2 ,m 2 ) 



(C6) 



and for p = this expression becomes 



£ m , m (°) = 7^9 m 6i(m 2 ) = G 2 (m 2 ). 
2m 



(C7) 
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